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Abstract 



Lemaitre-Tolman-Bondi models as specific spherically symmetric solutions of 
ly-j , general relativity simplify in their reduced form some of the mathematical ingredients 

of black hole or cosmological applications. The conditions imposed in addition to 
\& ■ spherical symmetry turn out to take a simple form at the kinematical level of loop 

quantum gravity, which allows a discussion of their implications at the quantum level. 
Moreover, the spherically symmetric setting of inhomogeneity illustrates several non- 
trivial properties of lattice refinements of discrete quantum gravity. Nevertheless, 
the situation at the dynamical level is quite non-trivial and thus provides insights to 
the anomaly problem. At an effective level, consistent versions of the dynamics are 
presented which implement the conditions together with the dynamical constraints 
of gravity in an anomaly-free manner. These are then used for analytical as well as 
numerical investigations of the fate of classical singularities, including non-spacelike 
ones, as they generically develop in these models. None of the corrections used 
here resolve those singularities by regular effective geometries. However, there are 
numerical indications that the collapse ends in a tamer shell-crossing singularity prior 
to the formation of central singularities for mass functions giving a regular conserved 
mass density. Moreover, we find quantum gravitational obstructions to the existence 
of exactly homogeneous solutions within this class of models. This indicates that 
homogeneous models must be seen in a wider context of inhomogeneous solutions 
and their reduction in order to provide reliable dynamical conclusions. 
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1 Introduction 



Black holes provide one of the most active areas of gravitational research, and a prime 
testing ground for quantum gravity. Many aspects can already be analyzed under the 
assumption of spherical symmetry, which is sufficient to describe non-rotating black holes. 
In vacuum, for instance, a reduced quantization is possible [TJ [2j [3] . This reduction has thus 
often been used in diverse approaches to quantum gravity, and also loop quantum gravity 
[U El E] has provided its own formulation [3, El El (101 ELI]- This allows one to use the loop 
representation constructed in the full theory in a simpler setting in which, as one hopes, 
one can find and analyze physical solutions. The result is a set of many coupled partial 
difference equations for a wave function on midisuperspace, which reflects the discreteness 
of spatial geometry realized in a loop quantization. As a consequence, spacelike classical 
singularities as they occur e.g. in the best-known case of the Schwarzschild solution are 
removed because their high curvature regimes no longer present boundaries to quantum 
evolution p2l[T3]. So far, these present the only inhomogeneous singularities analyzed by 
the methods of loop quantum gravity, and the situation appears much more non-trivial 
than in homogeneous models. Due to this complexity, several questions remain open. 
In particular, inhomogeneous models allow not only spacelike singularities but also null 
or timelike ones, which may even be naked; see e.g. [TH [TBI [TE] and references therein. 
Their fate presents additional problems of high interest which can be studied in spherical 
symmetry. 

Unfortunately, even the spherically symmetric equations in loop quantum gravity are 
difficult to tackle, and so it becomes interesting to look at further reductions which can 
preserve the physical setting but provide mathematical simplifications. Some of the possible 
effects have been introduced in a more ad-hoc manner in the models of [TT1 [T8l [T9l 120] , and 
were used mainly to understand implications on the horizon dynamics. A drawback of such 
an approach is that quantum corrections are not clearly linked to key ingredients of a loop 
quantization, or any quantization one could apply to the full, unrestricted theory. We thus 
intend to introduce new models which, starting from the spherically symmetric reduction 
along the general lines of [S] , provide further simplifications. Moreover, we will pay special 
attention to the anomaly problem which is required for full consistency. The possibility 
followed here is a reformulation of the spherically symmetric constrained system as it is used 
in the canonical definition of Lemaitre-Tolman-Bondi (LTB) models [211 1221 [23]. In this 
classical reduction, one solves the diffeomorphism constraint and, inserting the solutions 
into the Hamiltonian constraint, simplifies the set of equations. This eliminates some of 
the equations and removes several gauge issues. The remaining equations are then easier 
to quantize and solve. 

The imposed conditions make the constrained system partially second class, and thus 
are not straightforward to implement at the quantum level. Symmetry reduction, which 
is also partially second class and can be done at the kinematical quantum level [H [2U [253 
I2S1[2T], provides some guidelines, but the gauge fixing conditions required for LTB models 
are more complicated. While these issues are discussed here, we postpone a direct imple- 
mentation at the dynamical quantum level and rather perform in this article an analysis 
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of the quantum correction terms one can expect for the classical equations. Nevertheless, 
we will be able to highlight several important issues at the quantum level, such as the 
role of dynamical lattice refinements of discrete states j28j [29]. In particular, we will see 
that not all refinement schemes discussed so far in anisotropic homogeneous models can be 
embedded in spherically symmetric ones. But a large class still remains because only the 
direction dependence, not the size dependence of refinements is determined. 

What we provide in the main part of this paper is an analysis of consistent sets of 
quantum corrected equations of motion for a type of non-perturbative inhomogeneities, 
which has not yet been available in other models studied in the framework of loop quantum 
gravity. Previous equations either refer to homogeneous models [30], EH [32j [33] , where strict 
effective equations can be derived as analyzed in a canonical setting in [3U [35], [36] , but the 
consistency issue trivializes, or to perturbative inhomogeneities as in [371 ES EH HO, EH 
H2J H3j [H] . In our case, inhomogeneities can be non-perturbative but their equations are 
arrived at only after performing a partial gauge fixing. 

We will thus be able to restrict the form of possible quantum corrections by the con- 
dition that they provide consistent deformations of the classical equations. Thus, we will 
discuss how quantum corrections due to loop quantum gravity should occur in order to 
guarantee consistent LTB-like solutions. We will not be deriving strict effective equations, 
but we will determine consistent anomaly-free consequences of different types of quantum 
corrections. For several of the consistent choices we study consequences for inhomogeneous 
spherically symmetric systems based on analytical as well as numerical considerations. Our 
main interest here, as often in this context, is the fate of classical singularities. Compared 
to other models which by now have been studied extensively in loop quantum cosmology 
[15] , a new feature is the possibility of non-spacelike singularities in LTB systems. Interest- 
ingly, none of the corrections studied here resolve such singularities directly (nor spacelike 
ones in these models), in contrast to what their analogs do so easily in homogeneous mod- 
els. However, numerically we find indications that shell-crossing singularities generically 
occur in the presence of the corrections, where only shell-focusing singularities would do in 
classical theory. The shell-crossing singularities are expected to be weakly extendable, or 
to be avoidable by realistic matter. This is one example for new effects in inhomogeneous 
models, which not surprisingly prove themselves much more non-trivial than homogeneous 
ones. 

2 Spherically symmetric spacetimes and LTB vari- 
ables 

The line element in a general spherically symmetric spacetime with polar coordinates 
(x, <f) can be written as 

ds 2 = -e 2u dt 2 + e 2X dx 2 + R 2 (d$ 2 + sin 2 $d V 2 ), (1) 

where v = v(t,x), A = X(t,x) and R = R(t,x) are functions of the time and radial 
coordinates t and x. Alternatively, the form suitable for a canonical analysis, as done 
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below, is 



ds 2 = -N{x, t) 2 dt 2 + L{x, t) 2 (dx + N x {x, t)dt) 2 + R(x, t) 2 (d$ 2 + sin 2 fidip 2 ) (2) 

where N, N x , L and R are again free functions of the radial coordinate x and time t only. In 
a canonical quantization, the lapse function N and the shift vector N x appear as Lagrange 
multipliers of the Hamiltonian and diffeomorphism constraints, respectively. They are thus 
not dynamical, unlike the remaining functions L and R which have non-trivial momenta 
Pl and Pr, together forming the kinematical phase space of spherically symmetric gravity. 
In (pQ), the shift vector has already been chosen to vanish, while v and A are related to N 
and L in obvious ways. 



2.1 Matter source 



rpfj, 
U 



(3) 



We start with a spherically symmetric matter source whose stress-energy tensor is 
diagonalized in the form 

/ -e \ 
E 
n 

\ o o o n ) 

where e = e(t,x), S = £(£, x) and II = U(t, x) are the energy density, the radial pressure 
and the tangential pressure, respectively. Einstein's equation and the conservation law 
then reduce to a set of partial differential equations given by 



m = 4irGeR 2 R' 
m = -4nGi:R 2 R 
R! = Rv' + R'X 



E' 



( e + Ey-2(s-n)^ 



-2u 



R-vR 



,-2A, ,/ D ' 



i/R = - 



m 



-7r(E + 2IT) J R 



(4) 
(5) 
(6) 

(7) 



where we have introduced 



m 



2 



R ,2 e~ 2X 



R 2 e~ 2v 



(9) 



called the Misner-Sharp mass. 

The Misner-Sharp mass m is generally defined as 



m 



-R(l - V A RV A R) 



(10) 



where A takes values and 1 corresponding to the 2-manifold coordinatized by t and x. 
At spherical trapping horizons [16] we have R = 2m. (To see this, one can use the simple 
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criterion for the spherically symmetric trapping horizon as the place where constant area 
radius surfaces become null [57]. This is the boundary of the region of spherical trapped 
surfaces; there may be non-spherical trapped surfaces outside that region j48j [49].) In an 
asymptotically flat spacetime, the asymptotic value of the Misner-Sharp mass at spatial 
infinity equals the ADM mass Madm- We can define a conserved current from the Misner- 
Sharp mass, which is called the Kodama current and given by 

j A = e AB V B m, (11) 

where €ab is the antisymmetric tensor satisfying €ab£ B c = 9ac- By definition, the current 
j A satisfies the conservation law V a3 A — 0. The Kodama current provides the energy 
density 

m! 

(12) 



4ttGR 2 R' 
in agreement with 

2.2 Classical LTB spacetimes 

In classical Einstein gravity, there is an exact solution which describes a spherically sym- 
metric collapse system sourced by inhomogeneous dust: the Lemaitre-Tolman-Bondi solu- 
tion [2H [221 [23]. For dust, we have E = IT = such that Eq. (jgj) implies that 

1 

m = -F(x) (13) 

is an arbitrary function of the radial coordinate only. From Eq. (jlj), we then have 

F' 

C = 8ttGR 2 R> ^ 

for the dust density. Then, Eq. fl7|) implies that v = u(t) can be made to vanish by rescaling 
the time coordinate t such that it becomes the dust proper time. Eq. ([6]) can be integrated 
to give 

(.R'e^) 2 = 1 + K (x) , (15) 

where k(x) > — 1 is another arbitrary function. The resulting line element can be written 

as ^ 

ds 2 = -dt 2 + —, , dx 2 + R 2 (d$ 2 + sin 2 tfdip 2 ) (16) 
1 + k(x) 

where the only dynamical function left satisfies, from Eq. (J9j), 

R 2 = k{x) + . (17) 

Here we concentrate on the special case where k(x) = 0, called marginally bound LTB 
solution, such that 

ds 2 = -dt 2 + R' 2 dx 2 + R 2 (d$ 2 + sin 2 tfdip 2 ) (18) 
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and Eq. ([17)1 becomes 

R 2 = ^ • (19) 

At a trapping horizon, this expression equals one. We can integrate the equation of motion 
to 

2^3/2 

4 - t -M = ± 3FW*- (20) 

where t s (x) is an arbitrary function specifying initial values at t = 0. Its freedom can be 
absorbed by defining the function R in terms of the coordinate x at t = 0. For instance, 
requiring i?|t=o = £ and choosing the lower sign in ([2"U[) for a collapsing model, we have 
t s ( x ) = lx^/ 2 F{x)-^ 2 and thus 

/ Q , \ 2 /3 

R(x,t)= (x^- 6 -^F{x)t) . (21) 



2 v v , y 

(An example is the self-similar solution R = x(l — at/x) 2 ^ obtained for F{x) = \x as 
studied in [50].) Of special interest in those solutions is the small- x behavior of R which 
may be regular or singular. In particular, the energy density ( fbil) can be divergent where 
R' = 0, giving rise to shell-crossing singularities. While solutions are no longer valid beyond 
this point, shell-crossing singularities are deemed avoidable by more realistic matter models; 
see e.g. |51j . Moreover, solutions are extendable through shell-crossing singularities in a 
distributional sense [52j[53]. In addition to this, energy density as well as the Ricci scalar 
can diverge when R = 0, giving rise to a shell-focusing or central singularity which is the 
main interest here. The properties of shell-focusing singularities in LTB models have been 
extensively studied [5U [55J EH [57] . In particular, it was proven that naked singularities 
form generically in this system [55[ |58[ [59] . 

This provides a new arena for effects of quantum gravity, which we will come back to in 
more detail in Sec. [5j We will be using loop quantum gravity which requires Hamiltonian 
techniques. In the canonical formulation of spherically symmetric Einstein gravity, the 
Hamiltonians of the gravitational sector, H gTav , and of the matter sector, i^dust, for the 
LTB system provide the expressions 

R 2 R> + R(R 2 )' (R 2 R)' F 

-"grav — ~ 2G ' dust _ 2G 

after replacing the momentum in terms of R. The total Hamiltonian constraint -ffg rav + 
#dust = then yields Eq. ( fl9l) . Below, we will provide a detailed canonical analysis based 
on Ashtekar variables, which will allow us to incorporate some corrections as they are 
expected from loop quantum gravity. (See [601 ED E2] for a canonical analysis in ADM 
variables, and [63l IM] for an analysis in arbitrary dimensions.) 

2.3 Connection variables 

For a canonical formulation only the spatial part of the metric provides the dynamical 
degrees of freedom. Moreover, the most successful canonical quantization of gravity, loop 
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quantum gravity, is based on a densitized triad Ef rather than a spatial metric q a b, which 
are related by EfE% = &et(q cd )q ab . (We use tangent space indices a,b, . . . and internal gauge 
indices i, j, . . ..) The spatial metric for a spherically symmetric system in components of 
this variable is given by 



ds 2 = E *^A x 2 + \E x (x)\m 2 + sin 2 ^ 2 ) 

\E X [X '! 



(23) 



where instead of the spherically symmetric spatial metric components L (which equals 
E 9 j \/\E x \) and R (which equals yj^j) the spherically symmetric triad components E x 
and E v appear. Written as a densitized vector field taking values in su(2) with basis Tj, 
these components define the spherically symmetric densitized triad by 

E = E x (x)r 3 smi}— + (E 1 (x)r 1 + E 2 (x)t 2 ) sintf— + {E\x)t 2 - E 2 {x)r 1 ) — . (24) 
ox ov dip 

such that (E v ) 2 = (E 1 ) 2 + (E 2 ) 2 . (For more details on this decomposition see 0IB]- Notice 
that the sign of E x is not restricted to be positive, while E v is defined to be non-negative. 
The sign of E x thus determines the orientation of the triad since sgndet(i? J a ) = sgnE x ; it 
plays an important role for fundamental singularity resolution [13], but not for most of the 
analysis done in this article.) The remaining freedom of the three functions E x , E 1 and 
E 2 not contained in the components E x and E v is pure gauge since it does not affect the 
metric. It can be parameterized by an angle rj = &rcta.n(E 2 / E 1 ) which is subject to U(l)- 
gauge transformations. The corresponding Gauss constraint is G[ X] = J dx\(x)((E x )'+P r ') 
where P v is the momentum of rj. 

Fields canonically conjugate to the original triad components E 1 are Ashtekar connec- 
tion components of 

A = A x (x)r 3 dx + (Ai(x)ti + v4 2 (x)r 2 )d^ + (Ai(x)t 2 - A 2 (x)r 1 ) sin ^dy? + r 3 cos fidip . (25) 

Also here we introduce the U(l)-gauge invariant quantity A 2 V = A 2 +A 2 and the gauge angle 
(3 = SLYcta.n(A 2 / Ai) . Only the difference a := rj — (3 of the angles is gauge invariant, while 
each of them can be changed by the same amount with a gauge transformation. However, 
the new parameter A^ is not canonically conjugate to E v ; see for details. In terms 
of the U(l)-gauge invariant triad component E v , its conjugate variable is rather given by 
A^ cos a = jK v which happens to be proportional to an extrinsic curvature component [H] . 
(The constant of proportionality is given by the Barbero-Immirzi parameter 7 [65, 66J.) 
The canonical pairs we will be dealing with are thus (A x , E x ), {^K^, E^) for which 

{A x {x),E x (y)} = 2 1 G8(x,y) , {iK v (x),E*>(y)} = -yG6(x,y) (26) 

and another pair (rj, P v ) for the gauge angle. The relation between the remaining Ashtekar 
connection component A x and extrinsic curvature is A x = jK x — rj' where the spatial 
derivative v( of the angle is the (negative) x-component of the spin connection. Its angular 
components are given in terms of 

(E*l 
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which determines the relation between A 9 and as A 2 = V 2 + ^K 2 without reference 
to a. 

Solving the Gauss constraint removes the gauge angle and its momentum, and implies 
that invariant objects can depend on the x-component of the Ashtekar connection only 
through the extrinsic curvature quantity A x + rj = ~jK x . We are then left with only two 
canonical pairs (K X ,E X ) and (K^^E^). A further reduction can be implemented by using 
the LTB form of the variables corresponding to the metric ffl8l) . Comparing the spatial 
parts of the metrics ffl8l) and (1231) one obtains the LTB relation 

E«{x) = \\E x {x)\'. (28) 

In particular, we have T v = —sgnE x everywhere. For a complete reduction of canonical 
degrees of freedom, this is to be accompanied by a condition between K x and K v . We will 
determine this by consistency with the constraints, such that (1281) is preserved in time. 



2.4 Constraints 

Constraints in real Ashtekar variables as used here have been derived in [HI [9] . The diffeo- 
morphism constraint is given by 



D[N X ] = J dxN x (2A[E l + 2A' 2 E 2 - A X (E X 



(29) 



where A\ = A^sm/3, A 2 = —A^cos/3, E 1 = E^sinr], E 2 = —E^cosr/. Using these 
definitions in the above equation one obtains 



D[N X ) = J dxN x (2(A V cos a)' E* + 2(a' + (3')E V A V sin a - A X {E X )') 



(30) 



Using the relations A^cosa = jK v , A^siaa = T v , a + (5 — r), A x — —rj + jK x and 
r^, = — (E X )'/2E' P , this can be expressed solely in terms of U(l)-gauge invariant quantities 

as 

D[N X ] = J dxN x (2 1 K'^ - 1 K X {E X )') . (31) 

We can consider the marginally bound LTB condition, as given by (|28j) . as a gauge- 
fixing condition of the diffeomorphism constraint to form a second class pair. On the gauge 
fixing surface, the diffeomorphism constraint can be replaced by 

D[N X ] = J dxN x j{E x )'{sgn{E x )K' v - K x ) . (32) 

Thus, for the diffeomorphism constraint to be satisfied one must have 

K' v = K x sgnE x . (33) 
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This provides the LTB condition for the conjugate variables to E 1 analogous to (I28I) . 
Having solved the second class constraints resulting from gauge- fixing D[N X ], we proceed 
to the Hamiltonian constraint and insert the solutions. 

In spherical symmetry, the gravitational part of the Hamiltonian constraint is 

#grav[iV] = ~ J dxN(x)\E'T 1/2 ((1 - K + K%)E* + 2r x K v E%A x + rf) + 2E^) . 

(34) 

With r^, = —1 under the LTB condition in (I28I) . we have 

^ g rav[iV] =~J dxN(x)\E x \- 1/2 {Kl ELP + ^ K fE x K x ) (35) 
or, in unsmeared form, 

1 / K%E* . \ 

= i-^= + 2K v K x ^/\&\\ (36) 

where E^ and K v are to be understood as functionals of E x and K x , respectively, via 
the LTB conditions, or as functions of R and its momentum Pr. Thus, H gT&v is left as 
the sole constraint to restrict initial values and generate equations of motion for LTB 
models. In fact, the conditions (1281) and fl33|) can be seen to be consistent with the classical 
equations of motion for E x , E v , K x and K v generated by fl36l) (in addition to solving the 
diffeomorphism constraint identically). Thus, the reduction to LTB form is dynamically 
consistent and provides correct space-time solutions. Maintaining this consistency will be 
our main guideline in the analysis of quantum effects. 



3 Effects of a loop quantization 

Like any quantization of an interacting system, quantum gravity implies corrections to 
classical equations of motion and thus forces us to re-address consistency issues of classical 
reductions. As we will see explicitly, for a loop quantization this provides valuable feedback 
on the consistency or possible anomalies of the overall framework. In LTB systems, once 
a consistent reduction has been found, it can be used for applications to black hole singu- 
larity issues. We focus on loop specific issues which do not arise in Wheeler-DeWitt type 
quantizations of [60J ED E2] which have for instance been applied to Hawking radiation 

[BZ1ESIEI!. 

3.1 Basics of spherically symmetric loop quantum gravity 

A loop quantization [69] of spherically symmetric gravity is based on holonomies 

h e (A x ) = exp I ±ifA x dx ) , h v (K v ) = expfaK^v)) , ^(77) = exp(ir}(v)) (37) 
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for the configuration variables instead of linear expressions in connection or extrinsic cur- 
vature components. Here, we have used edges e and vertices v in the radial line. The 
use of these variables is strongly motivated by general results of loop quantum gravity 
[HEIE]: Background independence, i.e. quantization in the absence of a metric other than 
the physical one determined by Ef, is realized in a well-defined representation of smeared 
basic variables for holonomies together with fluxes as 2-dimensionally integrated densitized 
triads. But operators for connection or extrinsic curvature components themselves do not 
exist. 

An orthonormal basis of gauge invariant states in the connection representation is given 
by |B] 

T gAv = IT ex P ( fae I (Ac + 7]')dxj JJexp(z/i„7f^(t>)) (38) 

with integer labels k e and positive real labels fi v on edges e and vertices v, respectively, 
forming a finite graph g in the 1-dimensional radial line. The labels determine the connec- 
tion dependence by irreducible representations of the groups spanned by the holonomies. 
(These groups are U(l) for Ar-holonomies and the Bohr compactification Meohr of the real 
line for i^-holonomies; see [8] for details.) The densitized triad, i.e. momenta conjugate 
to the connection components, will be derivative operators which are quantized in the full 
theory in the form of fluxes f s E^r l n a d 2 y as integrals over surfaces S with co-normals n a . 
Their explicit action in spherical symmetry is 



i)1 ^e+(x) + k e -(x) , 



E*(x)T 9 , k! , = 7 4 w 2 W T 9 A, (39) 
/ E*Tg )ktli = 74y>^,^ ( 4 °) 

where £p = Gh is the Planck length squared and e ± (x) denote the neighboring edges to 
a point x, distinguished from each other using a given orientation of the radial line. (We 
have k e +M = fc e - (a) if % is not a vertex of the graph.) The -E^-operators only exist in 
smeared form after integrating over arbitrary radial intervals I. All flux operators have 
discrete spectra: eigenstates as seen in (1391) and (140p are normalizable. But only E x has 
a discrete set of eigenvalues, while .^-eigenvalues fill the real line. (Their eigenstates are 
elements of the non-separable Hilbert space of square integrable functions on the Bohr 
compactification of the real line.) 

These basic operators can be used for composite operators as well, providing well- 
defined but rather complicated constraint operators. Instead of dealing directly with these 
operators, we will extract some typical effects as phenomenological corrections to the clas- 
sical equations and analyze their potential implications with more ease. This has been done 
in quite some detail in cosmological applications [TQl [39l SQl SD EH S21 S31 SH [TH [721 [73] , 
and we start the same in this paper for inhomogeneous models in the spherically symmet- 
ric context. In this way, we provide the first examples where full inhomogeneities, rather 
than perturbative ones as in cosmology, are studied in this way. But we emphasize that 
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we do not derive complete effective equations as per [MJ [35J, [36] in this paper, which is 
rather an exploration of possible effects. Nevertheless, restrictions by consistency already 
provide interesting lessons. The general consistency issue is similar to that studied after 
partial gauge fixings in spherical symmetry in [TT], where it was analyzed based on the 
general ideas of [7U [75], and in [76]. Our analysis here provides complementary results 
in a different setting, where we make sure that potential anomalies arising from different 
correction terms cancel each other. Moreover, the effective treatment allows us to arrive 
more directly at properties of physical solutions. 

3.2 Implementing the LTB conditions 

The LTB conditions (|28|) and (I3"3l) in terms of densitized triad and extrinsic curvature 
variables are well-suited to an implementation at the level of spin networks. They refer 
directly to basic expressions of the quantization (provided that one just exponentiates the 
relation (1331) to result in holonomies) and can thus easily be formulated as conditions for 
kinematical states. In this way, the LTB reduction can be performed at the quantum level. 
However, consistency issues of the dynamics are not easy to deal with at the complete 
quantum level. We will therefore describe the kinematical constructions here, proceed to 
a phenomenological effective level in Sec. [4] and then study its consistency. By the link to 
the initial loop quantization, this indirect route will nevertheless provide feedback on the 
full theory. 

From the triad relation ([28]) we derive a condition for fluxes simply by integrating over 
arbitrary radial intervals X: 

^E*= l -\E x \ dx (41) 

where 81 is the boundary of X at which E x is evaluated, taking into account orientation 
to have the correct signs. This relation can be imposed on triad eigenstates (1351) . where 
([39D and p)} imply 

= ^(\k e +(v)\ ~ \K-(v)\) (42) 

for any vertex v. This directly eliminates all vertex labels in favor of the edge labels which 
remain free, analogously to the function \E X \ = R 2 which classically determines an LTB 
metric completely. 

On these reduced states, it turns out, the LTB condition for holonomy operators is 
already implemented. Upon integration and exponentiation, we have 

exp (^isgn(E x ) J(A X + r]')dxj = exp {\vyK v [v x )) exp (-^K v (v 2 )) (43) 

expressed solely in terms of elementary holonomy operators. This condition is realized in 
the sense that the left and right hand sides, as multiplication operators, have the same 
action on solutions to the LTB condition satisfying (j4"2"l) . In fact, the left hand side simply 
increases the label of the edge between v\ and v% (which we assume to be two adjacent 
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vertices) by one. Thus, it changes both fc e +(«i) an d k e -(v 2 ) by ±1 depending on their sign. 
The two operators on the right hand side, on the other hand, change the vertex label fx vi 
by | and fj, V2 by — ~ in the right way to respect the condition ( 1421) if it was realized for the 
original state. (If there are vertices v between v\ and t>2, k e +M and k e -r v \ change by the 
same value such that ()42l) remains implemented without changing fi v .) 

Notice that, unlike conditions for a symmetry reduction, the two LTB conditions 
for densitized triads and extrinsic curvature have vanishing Poisson brackets with each 
other (but not with the constraints). Thus, the curvature condition can indeed be imple- 
mented on the solution space of the triad condition. It does not add further conditions 
for states because they are written in a specific polarization. LTB states are thus sim- 
ply represented by a chain of integer labels k n for n = 0, 1, . . . which represents spatial 
discreteness (a 1-dimensional lattice of independent sites) as well as the discreteness of 
quantum geometry (integer k n as eigenvalues of the area radius squared). In a connec- 
tion representation, they can be written as T^(z ,Zi, . . .) = Y[ n z n n wnere the assignment 
n i— > z n := exp(|z f e ^K x dx) is a generalized LTB connection. 

While states can be reduced immediately to implement the LTB conditions, further 
conditions do result for composite operators because fj4"3]) must be used if the action of any 
operator is to be written on the LTB states where (j4"2l) has eliminated vertex labels. This 
provides reductions, e.g. of constraint operators, such that characteristic quantum gravity 
effects in loop operators can be carried over to constraints for an LTB model. 

3.3 Inverse triad effects 

The first effect we turn to arises from the required quantization of inverse powers of the 
densitized triad, such as E x in (1361) . There is no direct quantization of such an inverse 
since E x in a loop quantization (1391) has a discrete spectrum containing zero, and thus 
lacks an inverse operator. Nevertheless, one can use general techniques to arrive at a 
well-defined operator which reproduces (E x )~ 1 in a classical limit. This quantization is 
based on the Poisson relation between a holonomy and the volume, which one can identify 
with an expression for ^L-y . More precisely, we have 4 , K'yGsgn(E x )E ip / ^/\E X ~\ = {A x , V} 

for the combination of triad components appearing in the first term of fl36|) . where V = 
An J dx a/ I E x | E v is the classical expression for volume in spherically symmetric setting. In 
these expressions, we follow general constructions of the full theory [77]. 

When quantized, the connection component is expressed through a holonomy, and the 
Poisson bracket becomes a commutator. In order to stay as close to the full theory as 
possible, we use SU(2)-holonomies 

h x (A) = exp(r 3 / A x ) = cos(± / A x ) + 2r 3 sin(| / A x ) (44) 

which in their matrix elements provide the basic quantities (1371) . (Generators of SU(2) 
are Tj = —\i<Jj in terms of Pauli matrices Oj\ path ordering is not necessary for radial 
holonomies thanks to the symmetry reduction, which reduces the gauge group to an Abelian 
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one [7].) The commutator for these holonomies becomes 

h x [h~\V} = V- C os{\!A x )V C os{\!A x )-sin{\!A x )Vsm{\!A x ) (45) 
+2r 3 (cos(§ / A X )V sin(| / A x ) - sin(| / AJV cos(| / A 2 



and appears in the constraint in the form tx{r 3 h x [h x 1 , V]). This can be used in a quanti- 
zation of 

tr(r 3 h x {h-\ V}) = -tr(r 3 2 {/ A X ,V}) = ~J{A X ,V} ~ ^,{4^} 

where £o is the coordinate length of the edge used. When inserted in the Hamiltonian 
constraint, £q for all edges discretizes the integration measure dx. Eigenvalues can be 
computed easily from the basic action of holonomies and fluxes: for the operator 



EVSMEX) ~ l tr(r 3hx [K\V]) (46) 



'sin(i / A X )V cos(| / A x ) - cos(| / A x )Vsm{\ / A 



2n~fGh 



we have eigenvalues 




= 2v%>|/^| (y\J\k e +{v) + ^ e -(i,) + 1| - \J\K+(v) + K-{v) - 1|) (47) 

where f is the starting point of the interval X used as the edge in the holonomy. 

Looked at for all values of E x , eigenvalues of the resulting operator do not agree exactly 
with the classical function E v / ^J\E X \ but show deviations especially at small E x . We can 
parameterize this by a correction function a(E x ) as 

(48) 

where the subscript k(E x ) means that the eigenvalue of the operator is taken at label 
K+(v) + k e -(y) — 2E X /l^p as ^ follows from (I3H . Classically, i.e. for £ P — > 0, we have 
a(E x ) = 1, and this limit is approached by a for large E x . With the correction, the ex- 
pression for the inverse power of the triad component E x is finite, just as in the isotropic 
case [78j. The general behavior of the correction function is illustrated in Fig. [1] below. 
Similar constructions have been used before in spherical symmetry, see e.g. [17] . (While 
finiteness of inverse triad operators is realized in isotropic and spherically symmetric mod- 
els, this is not expected to be a general property [THl EDI- Nevertheless, inverse triad 
operators are well-defined in general situations of loop quantum gravity [81].) 
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One should note that E x refers to the total area of a whole orbit at radius x, which 
can have macroscopic values. In this case, a only slightly differs from one. However, 
a fully inhomogeneous quantization would refer to flux values of individual microscopic 
patches, where E x = ^2 n p n is a large sum of microscopic contributions p n , together giving 
the whole orbit area. The single fluxes p n are much smaller and closer to the Planck 
scale, which makes a differ from one more noticeably if these fundamental fluxes are 
used. This is an illustration of the fact that symmetric models often artificially suppress 
corrections from inverse triad operators, as first noted in [82]. Inverse triad corrections 
analogous to a have thus occasionally been underestimated, especially in isotropic models 
of macroscopic universes with large matter content. One can model the enhancements of 
corrections in fully inhomogeneous states even in symmetric models by using higher SU(2)- 
representations of holonomies in operators, not just the fundamental one as understood in 
(f46|) . Then, the expression for a changes essentially by replacing j£ P in (f48l) by j^ip if j 
is the spin of the representation. (See [531 154] for precise formulas in those cases.) For our 
qualitative analysis here we can focus on the expression (f4"8"|) . In fact, later applications 
mainly use the small- E x behavior near a center or a central singularity where corrections 
are strong for any j. 

Because classically, for £p — > 0, the function a(E x ) approaches one, the classical limit is 
correct if this function is inserted as a multiplier of l/\fE x in the Hamiltonian constraint. 
We can thus write fl5oD as 



= _ i ( •JW^\-JW ^5K KEf + KvK ^ 



(49) 



The form of corrections is usually not fully unique due to the presence of quantization 
ambiguities. Sometimes, however, they can be restricted more strongly by relating, when- 
ever possible, reduced expressions to what one expects in the full theory. This provides 
different options for specific corrections which can be analyzed for robustness and the phe- 
nomenology they imply. In the present case, the full theory may suggest an alternative 
corrected Hamiltonian 

H & = "2^7pf (i ^ + 2K * KxEX ^ ■ ( 5 °) 

Here, also the second term carries a correction function, which is motivated if one takes into 
account that the yf^J in the second term of fl36|) arises from a cancellation in E x j ^J\E X \ 
after inserting spherically symmetric variables into the full constraint. Thus, the second 
term could also be expected to have a correction by a. We will analyze both cases below 
and describe the differences they imply. 
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3.4 Holonomy effects and lattice refinements 

Another characteristic, and in fact eponymous, feature of loop quantum gravity is that not 
components of the connection but rather its holonomies are represented as operators. Since 
these are non-linear objects, additional corrections by higher order terms of the connection 
(or extrinsic curvature) will be present which again can be evaluated by including them 
as correction terms in phenomenological equations. This may appear as higher curvature 
corrections as they involve higher powers of extrinsic curvature, but we emphasize that 
this procedure will not give a complete picture since higher derivative terms are missing. 
These can be computed at an effective level [3H [35J, [36], which would require much more 
work not pursued here. 

We can correct for the holonomy effects in the Hamiltonian constraint (|36|) as follows. 
We assume K x to be fairly constant over a given edge (or part of an edge for graph-changing 
operators) of the graph whose holonomies appear as basic loop variables, so that JJ + K x « 
£qK x where £q is the coordinate length of the edge lying between the vertices t>_ and 
v + . Rather than using a precise loop quantization of ( 1361) and computing its expectation 
value in terms of holonomies, we make the following replacements in the Hamiltonian 
constraint: K v — > (-fS)^ 1 sxn^SK^) and K x — ► (7^o) _1 sin^-fC^o)- In addition to to, 6 
is a dimensionless parameter whose role is discussed below. With these corrections, the 
constraint becomes 



Hi" 1 ) = -—( si " 2 (7^) E * , sin( 7 ^) sin( 7 ^i ) r— - 
gmv "~ 1G\ 7 2 5 2 7 s 



(51) 



While this may not be the precise result from a complete effective constraint, it cap- 
tures the main effects of using holonomies as periodic, rather than linear, functions in the 
curvature components. Moreover, this simplest choice guarantees that also here the classi- 
cal limit, which involves a continuum limit 5 —>■ and i§ — > 0, is satisfied. Note, however, 
that one should not take the full functional form too seriously but rather view the sine 
functions as a convenient place-holder for a perturbative expansion in higher powers of the 
i^-components. Since there will be other corrections as mentioned above, they could easily 
dominate most of the expansion terms. 

Although the parameters S and £q appear in similar forms, their origins and roles are 
quite different. The parameter £ arises as the coordinate length of a radial edge along 
which we compute a holonomy. Its size is determined by the embedded graph we act 
on, as well as the precise form of the Hamiltonian constraint operator understood in the 
construction. While a Hamiltonian operator does not depend on the embedding and thus 
£o, that dependence would arise in the cause of computing an effective Hamiltonian as 
the expectation value in a state peaked on a classical geometry. Specifying the classical 
geometry requires one to partially fix the diffeomorphism gauge; the size of £q is then 
a direct measure for the discreteness of the state in this setting. Also 5 measures the 
discreteness, but it does not refer to the length of any edge. It is associated with curvature 
components along spherical symmetry orbits, and there is no room for orbital edges 
in this reduced model. To understand the meaning of 5, we again have to look at what 
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it should correspond to in a full, unreduced setting as we did in Sec. 13.31 to discuss the 
size of inverse triad corrections a. In a fully inhomogeneous setting, there would now be 
edges along spherical orbits whose lengths correspond to 8. For a configuration which is 
nearly spherical, there should be a regular distribution of edges forming a lattice on each 
symmetry orbit. The length of each edge, and thus 8, would decrease with an increasing 
number of lattice plaquettes J\f: 8 oc J\f~ 1 ^ 2 . In particular, for finer lattices we have 8^0 
just as £ ~~ > 0, approaching the continuum limit. 

The precise form of 8 depends on the exact state which is approximated by a spherically 
symmetric one. In particular, the argument shows that 8, unlike £q, may be phase-space 
dependent if geometrical growth is accompanied by a refinement of the lattice such that 
M{E X ) depends on E x , e.g. by a power-law form M{E X ) oc l^^. (Note that Air\E x (x)\ is 
the area of a sphere at radius x, and thus coordinate independent. In fact, a densitized 
triad E x in one dimension behaves like a scalar.) 

In this way, we are naturally led to a refinement scheme of phase-space dependent 
holonomies where point holonomies associated with depend on E x , while i^- holonomies 
along the inhomogeneous direction are triad independent. In a reduction to anisotropic 
but homogeneous models as in [82], this specific form has been shown to imply a dynamical 
law given by a fundamental difference equation which cannot be implemented equidistantly 
in minisuperspace variables. Equidistant versions of the difference equation, which would 
result if holonomies for a connection or extrinsic curvature component depend only on its 
conjugate triad component, i.e. E^ for i^-holonomies, are not embeddable in a spherically 
symmetric model. In this way, inhomogeneous models can reduce some of the freedom in- 
volved in choosing a refinement scheme for a homogeneous model. Nevertheless, at least 
a 1-parameter freedom of Af(E x ) oc |£ ,x | fe , or even a different functional behavior, is left. 
It is only the direction dependence of Af which is restricted, not the size-dependence. In 
particular, it is impossible to restrict the corresponding freedom in isotropic models. What 
we can also see is the fact that spherically symmetric refinement schemes which do not 
depend on any auxiliary scales can give rise to apparently scale dependent equations in 
a homogeneous reduction: when reduced to isotropy, a non-trivial refinement scheme can 
always be expressed by a function J\f(a) of the scale factor a. In contrast to E x , a is coor- 
dinate dependent since it rescales if spatial coordinates are multiplied by a constant. The 
reduction from spherical symmetry thus must automatically introduce a scaling dependent 
parameter f Q such that jV(a) depends only on the coordinate independent combination 
f a. This results in a well-defined way of non-trivial refinement schemes with all the free- 
dom as it is realized in spherical symmetry. Moreover, since the scaling dependence arises 
only in the reduction to homogeneity, it cannot be used as a reliable criterion to rule out 
refinement schemes if it is applied purely in homogeneous situations. 

The lattice refinement behavior is to be expected in any model on general grounds 
[28| [29] ; while a direct derivation of the behavior of M{E X ) from a full Hamiltonian op- 
erator is difficult, one can arrive at some properties and analyze their consequences phe- 
nomenologically. For cosmology, such work has been done in [7TJ [72J 1121 EI] and is initiated 
here for black hole physics. (The interior of the Schwarzschild black hole, which can be 
formulated as a homogeneous model, has been studied from this perspective in [82"ll8"5| l86].) 
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4 Corrected LTB models 



An LTB reduction at the dynamical quantum level of spherically symmetric systems is 
difficult because the combined algebra of LTB conditions and constraints, when seen as 
an extended constrained system, is not purely first class. Although, as mentioned, the 
classical LTB conditions are preserved by the equations of motion generated by the spher- 
ically symmetric Hamiltonian constraint, there is no simple off-shell algebra between these 
functionals which one could represent on the Hilbert space generated by spherically sym- 
metric spin network states. The various versions of quantum corrected Hamiltonians we 
have provided so far are thus not yet LTB reduced, although we have already removed 
the spin connection terms as they drop out in the classical LTB reduction. Moreover, 
although we did see that the classical LTB conditions can directly be taken over to the 
kinematical quantum level, such a step is much more complicated when combined with 
the quantum constraint algebra. If overall consistency with the constraints is required, the 
LTB conditions themselves may well require quantum corrections, too. In this section, we 
will be exploring the possibility of LTB-like solutions at a phenomenological effective level, 
allowing for corrections to constraints as well as the classical LTB conditions. 

At this stage, we still have two canonical pairs and two smeared constraints (the diffeo- 
morphism and corrected Hamiltonian constraint). The corrected constraints in this form 
are automatically consistent (i.e. first class) since the absence of the spin connection terms 
implies the absence of spatial derivatives in H grav ; the Hamiltonian constraints thus com- 
mute with themselves. Even if we add the non-dynamical dust contribution i^dust — F'/2G 
in terms of the mass function F(x), the system remains consistent. We have dropped the 
spin connection terms in anticipation of the imposition of LTB conditions, which solve 
the diffeomorphism constraint identically and thus show that the Abelian Poisson bracket 
{H[N], H[M]} = 0, realized even in the quantum corrected case, is correct. With such 
an algebra, the LTB form allows us to discuss the anomaly issue more easily. However, 
there is still a potential anomaly problem, whose solution allows us to draw feedback for 
quantizations of the Hamiltonian constraint: the reduction to constraints of LTB form 
is consistent only if there are LTB-like conditions, relating E v to (E x )' and K' to K x , 
which are preserved by the quantum corrected spherically symmetric constraints. The re- 
quirement of preserved LTB conditions will, as we will see, restrict the form of quantum 
corrections in different terms of the constraints. 

4.1 Consistent LTB reductions 

To derive equations of motion for the metric component R left as the only degree of 
freedom in an LTB metric, we eliminate K v and K x in favor of E v and E x using the 
equations of motion E 1 = {E^H} where I G {</?, x}. To evaluate these expressions we 
make use of the canonical Poisson bracket relations {K x (x), E x (y)} = 2G5(x — y) and 
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{Ktp(x), E v (y)} = G5(x — y) as they follow from (1261) . such that 

E x if- E X E V \ , s 

K v = — t= , K x = —== E* 52 

v 2^W\ VW\ \ 2EX J 

for the classical constraint. 

The key problem now is that (1281) combined with (1331) is no longer preserved by the 
evolution equations generated by fj49|) or (1501) for a ^ 1, as can directly be checked. If 
these equations were consistent, one could eliminate the variable E v in favor of E x in all 
equations of motion, which is then expressed as E x = R 2 . (From now on we assume E x > 
without loss of generality for the applications we are interested in.) Thus, a complete set 
of equations given by the Hamiltonian constraint for R and the second order evolution 
equation for R would be obtained. When the LTB conditions are not preserved, however, 
the Hamiltonian constraint equation for R will not be preserved by the evolution equation. 

Before deriving quantum corrected equations for R, we thus have to find LTB conditions 
suitable for the quantum corrected dynamics. The main conditions are (i) that they reduce 
to the classical LTB conditions when quantum corrections vanish, (ii) that they still solve 
the uncorrected diffeomorphism constraint identically as this is necessary for a consistent 
constraint algebra, and (iii) that they be preserved by the quantum corrected evolution 
equations. Condition (ii) is also motivated by the fact that finite diffeomorphisms are 
represented in loop quantum gravity directly by the action they generate on phase space 
functions without requiring any quantum corrections. 

4.1.1 Inverse triad corrections: first version 

In agreement with the diffeomorphism constraint, we make the ansatz 

(E x )' = 2f(E x )E^ , /\'. /(/-' )A r (53) 

where, as indicated, f{E x ) is assumed to depend only on E x in algebraic form. This func- 
tion will be determined by demanding that the new LTB conditions are preserved in time by 
( I49p . Writing the first constraint in the smeared form Cltb = / dxfi(x)(2f(E x )E lf — (E x y), 
we require that the Poisson bracket {Cltb, / dyH&lv}, which evaluates to 



vanishes for all n{x). Using fl53l) to remove the derivative terms we get 



(54) 



dz^(z) -UCyErVW^- + - * J ^Z . (55) 



df IfK^E* _ 2faK (p Ev 
dE x y/E^ \[E X 

For this to vanish for all //, the integrand must vanish which therefore gives the differential 
equation 

f{l-a) = 2E x 4l~ (56) 
M ' dE x v ; 
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for f(E x ) which, for a as in is solved by 



f(E* 



for E x > 7 4/2 and 



c l \fE x e~ (x l 2 



(v 7 ^ + y/E* - 74/2) (v 7 ^ + + 74/2 
Csv^exp (~\a + ~ arctan (y/E x /{i%/2 - W 

If + + 74/2 



1/2 



(57) 



(5f 



-> 00 the 
1 in the 



1 + 2~V4^' 



for E x < 74/2. The constant c\ is fixed by demanding that in the limit E 1 

corrected LTB conditions should go to their classical form. To ensure f{E x ) 

classical limit we have c\ = 2yfe. The functions / and a then have similar fall-off behaviors 

at large E x : f(E x ) ~ 1 + T 1 i%{E x )~ 2 + ■■■ while a(E x ] 

The second constant C2 is determined by continuity at E x 
2 5/4 e i/2-^/4 7 -i/4^-i/2_ It ig eagy tQ check that with thig form for ^ E x^ the other ltb 

condition, K' = f(E x )K x , is also preserved in time. 

We now eliminate the connection components from the Hamiltonian in favor of the 
triad components E x and E v using the equations of motion which give 



|/2, which gives c% 



E J 



2VE X 



E^ 



a- 



E X E* 
2E X 



(59) 



We then further eliminate E v using the new LTB conditions to obtain 



1 ( a(E x ) 2 (E x )' + E X (E X )' (E X ) 2 (E X )' 



2G \ 8/ (/•;' )•• - 2]\ r T?- 1. /•(/•:•' ) :i 2 

The total Hamiltonian constraint H^l v + H^ ust = then becomes 

a{E x ) 2 (E x y E x (E x y (E x ) 2 (E x y F , _ o 



8/(^)3/2 Af(E x f/ 2 



(60) 



(61) 



with the given ^-dependence of / and a. 

To obtain the evolution equation we take a time derivative of the first equation in (|59l) 
to obtain 

2E X E X - (E x ) 2 



K v 4|£*| 3 / 2 

Since K v is also determined by = {K v , HgH v } = —aK 2 j2\f& we have 



(62) 



4E X E X - (2- a)(E 



x\2 







(63) 
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1 



R/R* 



Figure 1: The correction functions a(R) (solid) and f(R) (dashed) where R is taken relative 
to R* := v^^p- 



We note that in this derivation of the evolution equation, f(E x ) does not appear anywhere. 
Nevertheless, its form is important for the mutual consistency of the evolution equation 
(|63|) and the Hamiltonian constraint equation (!6T|) . 

We can now write everything in terms of R using the relation E x = R 2 : 

= 2 VWTWM z VWEjm R (64) 

j£f> 

f 2R Cxp{ ^-«m) for J R> v ^72£ P 

— \ 2 5 / 4 R exp ( 1 (l-a(R))+ i arctan ( x /WJh& ,/2-R? j) -71-/4) . — ■ (65) 

' ^ VV P ^— -i for i? < v^72£ P 



7 1/4 v / 2FV«+V / - R2 +^p/ 2 

These functions are shown in Fig. [TJ The first order equation (I6TI) can be written as 

R 2 R'(a(R) - 1) + 2 j RAR' + i? 2 # = f(R)F' (66) 

and the evolution equation is 

2RR + R 2 + (a(R) - l)R 2 = . (67) 

These equations can explicitly be seen to be consistent upon using the differential relation 
( 1561) between / and a. 
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4.1.2 Inverse triad corrections: second version 



Starting from ( 1501) . which gives equations of motion 



E 1 



2aK v> VE* , 



aKu- 



E v 



, K 

n K oK 1 

V 2(E X ) 3 / 2 v 



/ E X 
da 



+ aK x VE' J 



E x dE 



Kt 



x H> 



E* 

s/E~* 



- 2 



da 



, , K V K X VE X 
dE x v 



-a- 



2v^ 



(68) 
(69) 
(70) 



we can proceed similarly. Also with these equations, the classical LTB conditions will not 
be preserved such that the reduction would be inconsistent. However, one can verify that 
the corrected LTB conditions 



£7* 



2a 



(E a 



aK x = K' 



(71) 



are preserved as before. 

Proceeding with these equations of motion and LTB conditions, we obtain 



R 2 R 



a.' 



F' = 



(72) 



as the Hamiltonian constraint equation in the presence of dust. Thus, the corrected equa- 
tion for R is 

R 2 R= a 2 (F(r) + c(t)). (73) 
The evolution equation derived via K v is 



2RR + R 2 = 2^R 2 
d log R 



(74) 



with a quantum correction on the right hand side. Taking a time derivative of (173]) . using 
F — and eliminating F from the resulting equation via (I73|) indeed produces (1741) provided 
c(t) = c = const. Thus, the system is consistent, and the freedom of c(t) in f[T3"j) is reduced 
to a constant which can be absorbed in the mass function. 



4.1.3 Holonomy corrections 

We now look for a consistent formulation with corrections due to holonomy effects. It turns 
out that the Hamiltonian (|5T!) leads to equations which are algebraically complicated to 
handle. To avoid technical difficulties, we first look at a Hamiltonian where K v appears 
via the function (7#) -1 sin('-f5K ip ) but K x has its classical appearance (i.e. the continuum 
limit £q — ► has been taken): 



rr(ma) J_ ^in 2 E* sin(« r~\ 

grav 2G\ 7 2 5 2 sf& !b xV J 



(75) 
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This describes holonomy corrections in regions where K v is large but K x remains small, 
or for states with a dense radial lattice such that £q is small. Again, the classical LTB 
conditions would not be preserved, and therefore we look for an alternative of the form 



(E*) , = 2g(K lp )E* , K'=g(K v )K x 



(76) 



compatible with the diffeomorphism constraint, where giK^,) is assumed to depend only 
on Kp. As before, this dependence will be self-consistently verified by demanding that the 
new LTB conditions are preserved in time: {J dx/j,(x)(2gE' p — (E x )'), J dyH^ly} = for 
all n(x). This Poisson bracket evaluates to 



dzfj,(z) 



sm 2 (j5K v ) E v dg sm(^SK v ) cos^SK^) E* 



2g- 



1 2 5 2 y/E* dK v 
- 2gcos^5K lf )K x VE* + 2cos(j5K (p )K' tp VE^ + 



sin( 7 «^) (E x )' 



7# 



/ E X 



(77) 



Substituting for (E x )' and K 1 from the corrected LTB conditions, we have 



dzfj,(z) 



sin 2 ijSKy) E^ dg 

2A2 



7 2 <5 



^diC 
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smdSKj cos(7^) E v 
7<5 \f& 
2gsin(j6K (p ) E* 



+ 



j 5 



'E x 



(75 



For this to be zero for all /x, the integrand must vanish which implies the differential 
equation 

sin (jSKp) dg 



dK, r 



2(cos(7^) - l)g 



solved by 



g(K v ) = c cos 4 (j6KJ2) 



(79) 



(80) 



The classical limit g — > 1 for 5 — > fixes the constant of integration c = 1. One can check 
that the other LTB condition is consistent with this choice for g{K v ). 

As before we now eliminate the connection components in favor of the triad components. 
The equations of motion give 



E x 
E v 



E 2 



2 sin (•ySK lf 
7$ 

sin {^8K 9 ) cos {ySK v ) E v 



l5 sJe? 
We use (j8"Tj) to express the (co)sine function in terms of E x : 



+ cos {^K^) K X V& . 



(81) 
f82l 



sin (j5K„ 



7<5 E a 



/ E X 



cos ('ySK if 



±11- 



7 2 5 2 (E X Y 



(83) 
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and solve (1521 for 



sin (7^) cos (7^) 



K x = = + ^ v ' - w . (84) 

After substituting for the sine and the cosine (choosing the plus sign in the cosine), this 
can be written as 



Km =[&- - ^JlI ) (e* _ ffm'Y' 2 (85) 



2E X V 4 E 2 

Substituting for K x and sin (jSK^) back in the expression for the Hamiltonian ( 1751) we 
have 



H (ma) = _ J_ f E*E* _ {E*)*E* 

2G \Je* - h 2 5 2 (E x ) 2 4 ^) 3/2 



grav 



i6) 



The new LTB condition £ v = (E x )' /2g = (E x )' '/2 cos 4 (75^/2) allows us to eliminate £ ¥ 
using 



cos 4 ( 7 ^/2) =l Ul- I (87) 



which gives 

2(£ 3; ) / 



2 ' 



4 ' " 

Substituting for E 9 and its time derivative in (|86|) implies 

2E X (E X )' (E X ) 2 (E X )' 



(1 + - h 2 * 2 ^) V 1 -^ 2 ^ 2(^)i (i + v /i-h 2 ^)' 

1 2 b 2 E x [E x )\2E x E x E x - (E x ) 3 ) 

2(^)1 (l + ^1-I 7 2 ^) 3 (l - h 2 ^ 2 ^) 

We now derive the evolution equation consistent with fl86l) . On the one hand, we have 
the equation of motion 

K =-\K H (nia) } = s j^MM 1 (g9) 

W^grav / 2 ^ 2 ^2 V "^ 

and, on the other hand, differentiating ( 1511 with respect to time gives 



K= l - ( -jl- ^ 1 (90) 

* 2(308(7^)^^ 2(E X ) 3 / 2 ' 
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Combining these two equations and writing everything in terms of R, we have 



2RR + Ryi - 1 2 5 2 R 2 = 0. (91) 

The evolution equation can now be used to eliminate the second time derivative of 
E x from Hg V H a \ which together with E x = R 2 and in combination with the matter part 
provides the Hamiltonian constraint equation 



4R 2 R'Jl - j 2 5 2 R 2 + 8RRR' = F' 1 + J 1 - 1 2 5 2 R 2 J 1 - 1 2 5 2 R 2 . (92) 



We note that in the limit 6 — > we recover the classical equation which also justifies the 
choice of plus sign in ( 1531) . 

Finally, we could use a Hamiltonian where only K x has been replaced by periodic 
functions, 



BgP = ~± ^ + 2K^^^VeA . (93) 



1 ( K l EH> , « v sin (7^ o; 
tgrav 2G V " 74 

In this case, however, the corrected LTB conditions will take a more complicated form 
because correction functions will have to depend on all the phase space variables, as one 
can check by making an ansatz as before. We leave this complicated case open for future 
work and proceed with a general discussion and applications of the consistent versions 
found. 



4.2 Discussion 

We have provided several cases of consistent equations of motion for the variables of a 
metric of LTB form, but with dynamics carrying corrections as they are expected from 
loop quantum gravity. While we have discussed inverse triad and one form of holonomy 
corrections separately, they can be seen to be combined consistently simply in a multi- 
plicative form of the correction functions in the LTB conditions. For instance, the first 
version of inverse triad corrections and the holonomy correction we used can be consistently 
combined in this way to result in a Hamiltonian constraint equation 

aR 2 R ' + 2RM ' = f/jF > (94) 
yjl --i 2 5 2 R 2 

where f$[R] = f(R) f 1 + \j 1 — 7 2 <5 2 i? 2 ^ together with the evolution equation 



2RR = -aR 2 V 1 - -f 2 5 2 R 2 . (95) 
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With the second version of inverse triad corrections, we have 



Aa 2 R 2 R'\ll - - Aa 2 R 2 R' - ^ 2 S 2 R 4 R' + 8a 2 RRR' I 1 



a 2 



8a 2 R 2 R 2 R' / 1 2 5 2 R 2 \ _ 4 / 1 2 5 2 R 2 

, I 1 + 1/1-^2—1 =F'« 4 I 1 + 1/1 - 



& - ( 7 ^/2) 2 





and 



2ffli = 2J? _*^_4_^. (97) 
Vfl 1 - (74/2) 2 V " 2 

While general properties of an LTB reduction allowed us to keep the constraints consis- 
tent without severe limitations on quantum correction functions, consistency conditions did 
remain. The remaining constraints automatically form a first class system provided that 
the Hamiltonian constraint is free of spatial derivatives, which is realized if spin connection 
terms (or the 3-dimensional Ricci curvature) drop out as it happens under the classical 
LTB conditions. The consistency conditions arose at the level of formulating the LTB con- 
ditions, because the classical ones are no longer preserved under evolution corresponding 
to quantum corrected constraints. We thus corrected the LTB conditions, too, such that 
in their new form they were preserved under the quantum corrected equations of motion 
as they are generated by a Hamiltonian whose spin connection contribution vanishes. 

With these conditions we are still identically satisfying the classical diffeomorphism 
constraint and thus no new anomalies in the constraint algebra can arise. However, the 
corrected LTB conditions do not make the classical spin connection T v equal —1, which was 
assumed in the simplified classical Hamiltonians such as (149)) . Thus, to be fully consistent 
we must assume that the expressions containing the spin connection themselves carry 
quantum corrections and read f(E x )T^ = —(E x )'/2E ip for the equations in Sec. 14.1.11 
a{E x )r { " ] = -{E X )'/2E* for Sec. EJ and g(K <p )r$ IIa) = -(E X )'/2E* in Sec. I3T31 
Thus, for consistency additional corrections of this form must arise in the terms of the 
Hamiltonian constraint containing the spin connection in such a way that they vanish 
after imposing the quantum corrected LTB conditions. 

That the spin connection terms carry their own corrections is a reasonable expectation: 
There are inverse triad components and, in a loop quantization, the spin connection is 
rather indirectly expressed via A\ and = A l a — T\ oc {A l a , {H^, V}} using the Eu- 
clidean part of the Hamiltonian constraint [77]. Corrections are thus expected from 
the inverse triad as well as from holonomies. The specific form is difficult to determine 
because the full theory does not provide operators for the non-covariant spin connection 
components, but as demonstrated here it can be derived and justified by the production 
of a consistent set of equations. In fact, if such corrections occur, our equations provide 
a fully consistent LTB system. In this way, consistency determines what further quantum 
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corrections must be entailed by a primary correction such as a. Since not all corrections 
in a Hamiltonian can equally easily be computed, independent consistency considerations 
provide useful relations between different terms. For instance, the spin connection is more 
difficult to quantize than l/E x , and its corrections can thus more easily be found via 
consistency. 

An important physical implication is that this suggests additional effects because the 
space-time metric ( 118)) is no longer just corrected by different solutions for R(t) solving 
the corrected constraint and evolution equations, but also by an additional pre-factor in 
terms of a, f or g in front of L 2 oc (-R') 2 which is no longer exactly {R') 2 - This would, for 
instance, affect the appearance of horizons. 

5 Applications 

Our focus in this paper for applications of the above equations is the fate of the classical 
singularity which appears at R = 0. There are two possibilities for how such a singularity 
could be avoided in effective equations. Dynamically, R(t, x) may be bounded away from 
zero for all x, in which case the behavior shown would be comparable to a cosmological 
bounce. This can sometimes occur if isotropic cosmological models exhibiting a bounce are 
matched to a spherically symmetric outside region in a generalized Oppenheimer-Snyder 
manner [87]. The outright spherically symmetric situation studied here is, however, subject 
to different corrected equations and so one has to provide a new analysis. 

The second possibility is that the value R = is assumed, but that this does not result 
in a singular space-time just as Minkowski space in polar coordinates has R = at x = 0. 
If R = occurs, one thus has to proceed with a more detailed analysis to understand the 
space-time neighborhood of the region where R = 0. 

Compared to homogeneous equations, this problem is of a new quality. As we have 
seen, there is a non-trivial anomaly problem which we were able to resolve in different 
versions of quantum corrected LTB models. The presence of consistency conditions, which 
do not arise in homogeneous models because they are subject to just a single constraint, 
makes the form of quantum corrections more restricted. Thus, several different terms in 
the constraints must receive quantum corrections in a way closely related to each other. 
Still, we have explicitly shown that non-trivial quantum corrections are allowed. 

In addition to the anomaly issue, spatial inhomogeneity allows different types of sin- 
gularities in classical general relativity. In particular, not just spacelike singularities can 
occur as in homogeneous models, but also null [57)1 158] or timelike ones [SB] • This has in- 
teresting general ramifications concerning the consistency of quantum gravity in the sense 
of allowing stable ground states, as discussed in [89], and it underlines the interest in in- 
homogeneous models. In what follows, we present an initial analysis based on analytical 
as well as numerical methods. 
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5.1 Analytical properties 

If there is a "bounce" where the area radius R attains a non-zero minimum value, we 
have R = which can be substituted in the above equations to check the possibility for 
this to happen in quantum gravity. From Eq. ( 1661) (for the first version of inverse triad 
corrections) or Eq. fl92l) (for holonomy corrections) as well as the two combinations fl94|) 
and fl96l) we can immediately see that this is not possible unless we drop the condition 
F' > which classically avoids shell-crossing singularities. Thus, we either have to drop 
this condition, possibly taking into account quantum geometry corrections in the matter 
sector, or retain the non-bouncing behavior of the classical models. For the second version 
of inverse triad corrections, Eq. (!73|) . we would require F = at the bounce, which looks 
difficult to achieve in a generic collapse model. (Conditions on F may be avoided if R' 
diverges where R = 0, but this does not appear generic.) 

There does not appear to be a simple conclusion about bounces as they occur, e.g., 
in homogeneous models. We are looking at specific regimes and certain types of quantum 
corrections which, by themselves, may make a bounce difficult to occur. Moreover, we 
have restricted the analysis to marginal models, which classically includes spatially flat 
Friedmann-Robertson-Walker models as the interior region of Oppenheimer-Snyder col- 
lapse, but not isotropic models with positive spatial curvature. The latter (or scalar matter 
with negative potential |90j) would be required for a bounce based on inverse triad correc- 
tions [91J. There is thus no contradiction with known matching results based on isotropic 
interiors [87] , but the fact that a bounce does not follow straightforwardly, compared to 
the relative ease by which this can be obtained in isotropic models, may also be taken as 
a warning sign concerning the robustness of homogeneous bounces. 

Similarly, we have considered holonomy corrections only of a special form which made 
the analysis more manageable. Holonomy corrections give rise to bounces more generally 
than inverse triad corrections [921 [931 ESS E2] . One could thus expect that a full treatment 
of holonomy corrections should give rise to general bounces also in LTB models. However, 
even though we did not do such an analysis, such a bounce cannot be generic for the 
following reason: simply choosing a fine spatial graph and thus small enough £ Q makes the 
holonomy corrections studied here the relevant ones. Since these corrections do not provide 
an automatic bounce, a bounce cannot be generic in this inhomogeneous system. Finally, 
there is a third effect due to quantum variables which provides corrections in effective 
equations. Also this has not been included here, but it is unlikely to result in a general 
bounce given that it does not do so in isotropic models (where it could even prevent a 
bounce which would otherwise occur based on holonomy corrections [521 [55]). We thus 
conclude that singularities in LTB systems do not appear to be resolved by bounces. 

5.1.1 Corrected LTB equations as cosmological models 

As the simplest case, we first consider a vacuum solution where, in the absence of dust, 
F' = must be satisfied. If R' ^ holds, R = R(x), i.e. a static configuration, is a trivial 
solution to any of the corrected equations (RjoT) . f[75|) and f[9"2"j) . While this corresponds 
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to the classical Minkowski space solution, since R(x) can then easily be introduced as a 
coordinate instead of x, there are quantum corrections for small R: Our corrected LTB 
metric, using L = R' / f(R) for (l66l) to be specific, reads 



ds 2 = -dt 2 + 



dR 2 
f(R) 2 



+ R 2 (d$ 2 + sin 2 $dip 2 ) 



(98) 



which asymptotically presents Minkowski space. As we will discuss in more detail below, 
the appearance of f(R) shows that quantum effects originating in the spatial discreteness 
of loop quantum gravity spoil some of the exact symmetries such as spatial homogeneity 
known to exist in classical solutions. 

If there is dust, it is of interest to see whether we can have a Friedmann solution in 
this system. For this we choose x such that it coincides with the circumferential radius 
at t — 0, and make an ansatz of the form R(t,x) = a(t)x. If such a solution exists, as it 
does in the classical case, the LTB metric reduces to a Friedmann-Robertson-Walker one 
where a is identified with the scale factor. The dust density profile then becomes 



where €q is the initial uniform density. (According to our general choice of R(0, x) = x, 
the scale factor is normalized to a$ = 1 at t = in the cosmo logical context.) Substituting 
this into Eq. fl66|) . we have 



The left-hand side depends only on t while the right-hand side depends non-trivially on 
ax. Hence, the corrected LTB system does not admit a solution of Friedmann form. 

There is an additional effect which prevents Friedmann solutions for the corrected 
equations, because our LTB form of the metric receives quantum corrections, too, as a 
consequence of consistency. The corrected LTB metrics have coefficient L = R'/f(R) in 
the case of and L = R' /a(R) in the case of H^ n \ This changes the metric in addition 
to the corrected dynamics of the metric component R. In particular, the metrics are no 
longer homogeneous because of the non-trivial independence. If we were interested in 
spatial volumes of finite regions in constant t slices, for an approximate solution of the 
form R = a(t)x they would become V = 47ca 3 J dxx 2 / f(ax) and V = 4na 3 J dxx 2 /a(ax), 
respectively. This illustrates an interesting difference between these two cases which both 
come from inverse triad corrections: for f(ax), we have the small-x expansion f(ax) oc 
ax+0(a 2 x 2 ), while for a it reads a(ax) oc a 3 x 3 +0(a 4 x 4 ). Thus, in the first case the spatial 
volumes vanish at a = as in the classical case, while the second case implies diverging 
volumes V ~ J x~ 1 dx even of finite regions near a = 0. This suggests further implications 
of the behavior near a classical singularity, which due to the required inhomogeneity do 
not appear easy to discern. 




(99) 




(100) 
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5.1.2 Effective densities 



As a further consequence of corrections, we note that the mass function F(x) is no longer 
directly related to the Misner-Sharp mass. We need to distinguish the latter from the dust 
mass which can be defined as 

Fix) 

M{x) = -±L. (101) 

The asymptotic value Md us t = linx^oo M(x) corresponds to the total mass of dust. The 
Misner-Sharp mass, on the other hand, now takes the form 




in - ( 1 ( — ) + R- j (102) 

whose expression changes for the corrected LTB conditions because this affects the relation 
between R' and L. From the corrected mass, we can then derive an effective density 
e = m! / AnGR 2 R' . For instance, for the condition consistent with H&l y we have a Misner- 
Sharp mass 

m W = 1^(1 - f + R 2 ) = m class - ±R(f 2 - 1) (103) 
which, upon using fl66l) . leads to an effective density 

= (^§^ - HR) - l)(R 2 - 2f{Rf) - {f{Rf - 1)) . (104) 



Similarly, for the equations following from HgUl, we have m = \R(1 — a 2 + R 2 ) and 
thus 



Mi) 1 a 2 F> a 2 -l 2daR 2 a da \ 

~8nG\R 2 R' R 2 adRR RdR)' 1 ; 

In particular, in this case the horizon condition 2m = R reads R 2 = a 2 , which by the 
Hamiltonian constraint equation agrees with a 2 F/R. Thus, in terms of F the horizon 
condition F = R is uncorrected in this case, although R as a function of time is corrected 
compared to the classical behavior. 

The correction terms to the effective densities may be nonzero even in vacuum regions 
devoid of dust. Depending on the regime, they can be positive or negative according to 
the signs of a(R) — 1, f(R) — 1 and their derivatives involved. 



5.1.3 Existence of self-similar solutions? 

The classical equation can rather easily be analyzed using self-similar solutions; see e.g. 
|50j . One can first write the classical constraint equation as R 2 = F(x)/R and then, for the 
special case of a linear mass function F(x) = Xx, find an explicit solution for R(x, t) of self- 
similar form which depends on t only via the function 1 — at/x with a constant a = §\/X 
For such a self-similar solution, the structure of the singularity has been analyzed in |94j . 
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If we use the second version of inverse triad corrections, this equation is simply changed 
by multiplying the mass function with a(R) 2 . Thus, for a linear mass function there is no 
longer a self-similar solution. One would have to incorporate the new factor by changing 
the mass function, if a self-similar solution is to be obtained. But this is not straightforward 
since a depends not on x but on the unknown function R which is to be solved for. 

For the first version of inverse triad corrections the equation changes more radically. 
In this case, we can bring the constraint equation to the form of an integro-differential 
equation for R, 

(f ) 2 = 5 / d *nR(z))F'(z) - i/dfl(a - 1) (f )" . (106) 

If a — 1 is small, one can solve this iteratively by inserting the equation for (dR/dt) 2 in 
the integral: 

(^) 2 = ^ + | / Mf - l)F'(x) - i J dR(a - l)F(x)/R +■■■. 

The difficulty in solving this is that F(x) depends on x rather than R, so we have to know 
R(x) as a solution and invert it before doing the integration. But this equation already 
shows qualitatively that the quantum correction makes the solution more non-local, which 
may prevent the existence of self-similar solutions. Moreover, there will be additional time- 
dependent effects which do not occur classically. This is so because R(t, x), which we need 
to know in order to replace x by R in the integrand, also depends on t. Thus, the integrals 
are really time dependent, which one can understand as replacing the classical F(x) by a 
new function 

n*, t) = ^t/^-w./^-wp ^ dor, 

- J dR(a - l)/R J dx(f - l)F'(x) + j dR(a - l)/R J dR(a - l)F(x(R, t))/R +■■■ 

appearing on the right hand side of (dR/dt) 2 = J r (x,t)/R. 

This refers only to the case where a — 1 is small, i.e. we have perturbative correc- 
tions to the inverse triad effects. It would not allow one to analyze the deeper quantum 
regime where a differs significantly from one. For this regime we would have to use other 
techniques, such as the expansions of the following subsection. 

However, generally speaking, since the characteritic length scale £p is explicitly in- 
troduced into the corrections, we cannot expect this kind of self-similar solutions, which 
are called complete self-similar solutions or self-similar solutions of the first kind. With 
the characteristic length scale, we can only expect incomplete self-similar solutions, e.g. 
kinematic self-similar solutions in this context; see [95 | |96| [97]. 
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5.1.4 Small- a; expansion 

For R 3> y/^ylp, a — > 1 and the classical limit is recovered. However, for R < ^/yip, 
the deviation from classical theory becomes of order unity. Here, deep quantum effects 
might be revealed by a closer analysis. We point out that such a deep quantum regime is 
less reliable if only one type of quantum effect is considered. Nevertheless, an analysis of 
single effects can provide various possibilities and guide further developments. Moreover, 
the corrections studied here, based on inverse triad and holonomy corrections, can be 
combined without changing the conclusions. 

To have a regular center in an inhomogeneous case we assume that F and R admit the 
following expansions at the center. 

F(x) = F 3 x 3 + F A x 4 + ■■■ , R(t, x) = R!(t)x + R 2 {t)x 2 + ■■■ (108) 

where the dots denote higher order terms with respect to x, and Fj are constants but 
Ri may be t-dependent. In this way, the classical expression for energy density gives an 
expansion of the form 

e(t, x) = e (t) + e x {t)x + e 2 {t)x 2 + ■■■ (109) 

from F(x) and R(t,x), where e = 3F 3 /87rGi?f. Classically, the lowest order then gives 
Rl = F$/Ri with solution 

3^ n 2/3 



R 1= ^C±^tj , (110) 

where C is an arbitrary constant. Hence, for the collapsing case, R\ monotonically de- 
creases and becomes zero in a finite proper time — the central singularity develops where 
e — > oo. When we choose the radial coordinate x so that R = x at t = t , we find 

Rl = l±%-A(t-t )\ , (111) 



and hence, R\ vanishes at t — t s , where 

t s = toT^=- (112) 

This behavior can be checked also in the presence of quantum corrections to see if anything 

of the singularity changes. For the first version of inverse triad corrections, we use the 
small-x behavior 

/ O \ 3/2 /cpl-TT/2 

a =(-^) R l xS > f=\ H2— R ^ x ( 113 ) 



of the correction functions. As a result, if F(x) has a cubic term F^x 3 as the lowest order, 
we find i?i = 0. However, the additional factor of x in the dust energy density proportional 
to f(R)F'/R 2 R' then shows that we can allow a quadratic term F 2 x 2 in the mass function 
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to produce the desired regular expansion for energy density, although the total effective 
energy density is still diverging at the center because of the "vacuum" contributions in 

(HMD. 

Using the various series expansions in (1661) . to lowest order in x we get 



which is solved by 



Rl 



Rtit) = 1 ± 



2 3 / 2 e^-fF 2 
(74) 1/2 



2 3/2 



(74) 1/2 



1/2 



(t ~ t) 



(114) 



115) 



where the plus sign corresponds to a collapsing dust cloud and where we have chosen the 
initial condition Ri(to) = 1. From here we see that the central singularity, corresponding 
to Ri(t) = is formed at 



t = tn 



-i -1 



(74) 1/2 



116) 



For the sake of comparison we note that classically the central singularity forms at t = 
t + 2/3y/F^. In terms of the initial density profile, F 3 = 87rGe (0)/3 (note that e (0) is not 
the complete dust density profile at the initial time but a coefficient in the series expansion 
for the dust density) and therefore in terms of the initial density the time for singularity 
formation is t — to + 1/ ^/QixGeq. For the quantum corrected case a similar consideration 
gives F 2 = 47rGeo(7^p) 1//2 /2 3//2 e2~^ which implies that the central singularity forms at 
t = t + l/^f^Ge a . 

For the second version of inverse triad corrections, the equation to be solved is R 2 R = 
a 2 F from (I731 . For small R, (11131) implies that the lowest order term on the right hand 
side goes as x 9 if we use the same form of the expansion of F(x) as classically. In this case, 
we can show R\ = as in the first version. However, now our density is e oc (a 2 F')/R 2 R' 
if the divergent "vacuum" contribution is subtracted, which means that o?F should be 
required to have a leading term cubic in x, such that F itself can have lower order terms. 
In this case, the singularity is not prevented either. Nevertheless, a finite neighborhood 
may look different from the classical space-time near a central singularity which, however, 
would require a more detailed analysis. 

For our version of holonomy corrections, the equation to be solved is (JS2D - After 
expanding in powers of x and equating the coefficients on the two sides we obtain (to order 
R\R\ = F3 which is the same as in the classical case and gives 



x 



Ri(t) 



-VFsit-to] 



2/3 
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implying that the singularity occurs at time t = to + 2/3v^3 (the same as in the classical 
case). At order x 3 we have the equation for R2 which also does not have any quantum 
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corrections. Effects of the 7 2 <5 2 factor occur only at order x 4 and higher. Thus the small x 
behaviour with this correction is the same as for the classical LTB model. The combinations 
of inverse triad and holonomy corrections in ( 1941) and ( l96l) lead to the same conclusions. 

We conclude that there is no indication that the corrections implemented here prevent 
the LTB singularity from forming. In particular, naked singularities as they appear in 
these models do not seem resolved automatically by a loop quantization. Whether they 
are indeed naked singularities in the presence of quantum effects requires further analysis 
of the effective space-time: the surroundings of the singularity may be sufficiently different 
from the classical naked case such that the singularity becomes spacelike. However, we have 
shown that a situation which gives rise to a naked singularity classically also gives rise to a 
singularity (of some form) under the quantum effects considered here. An analysis whether 
naked singularities remain naked may be of interest in the context of cosmic censorship, 
which we will come back to in the numerical analysis. 

A correction not considered here is the effect of i^-holonomies which are computation- 
ally more complicated and may be crucial in some regimes. In fact, taking the general 
form (T2T1) of classical solutions indicates that K v = R = \JF/R is subdominant to 



near x = 0. This may present a chance for holonomy corrections to remove the singularity, 
after all. However, the size of holonomy corrections is also state-dependent: a small £0 sup- 
presses ^-corrections even if K x is large. Thus, as already noted these corrections cannot 
result in generic avoidance of singularities. (Note that also in [76J only ^-corrections 
were included, although the conclusion was that the Schwarzschild singularity could be re- 
solved in this way.) This clearly shows the non-trivial behavior of inhomogeneous situations 
compared to homogeneous ones. 

5.2 Numerical analysis 

Although some analytical results are available, the region of intermediate values for R 
not in the asymptotic regimes R -C a/t^p or R 3> a/7^p is difficult to analyze. We thus 
complement the preceding analysis by numerical studies of inverse triad corrections. 
For the first version, we transform the constraint equation with 






to 



K' 





(120) 
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We spatially integrate the constraint equation at each time step and evolve R using (11201) . 
The constraint is always satisfied within 10~ 4 accuracy. The spatial integration is done 
using a fourth-order scheme while the time evolution uses a second-order scheme. 

The second version is given by a set of equations which is first integrated by RR 2 = a 2 F. 
The numerical implementation of this equation is easy because we no longer need spatial 
integration. Our numerical scheme is second-order, which is sufficiently accurate and stable 
for the present purpose. 

It should be noted that there is scale invariance in both sets of equations which are 
invariant under scaling R(t,x), F(x) and t as R(t/j3,x), ft 2 F(x) and t//3, where ft is a 
positive constant. This scale invariance greatly simplifies the analysis. If the functional 
form of F(x) is the same up to a constant factor, the evolution is similar up to the scaling 
of time. We thus do not need to investigate the full parameter space. Through this scaling, 
the classical Misner- Sharp mass m c i ass scales as ft 2 m c \ ass , while the corrected mass m does 
not. For holonomy corrections, the scaling behavior is lost. Here, a dedicated analysis of 
the whole parameter freedom is required to draw reliable conclusions, which we postpone 
to future work. 



5.2.1 Initial condition 

The function F(x) corresponds to the conserved mass because it is constant for a comoving 
observer. Then, 

F' 

Pcons = 8nGR 2 R> (121) 
is regarded as the conserved mass density. Note that this differs from effective energy 
densities which incorporate quantum effects. 
We choose the radial coordinate x so that 

-fi-Ol^maxJ 

i.e. < x < 1 for the region of computation, where Rq(x) = -R(0, x). As matter models, it 
is useful to consider two different cases: (i) uniform models where 

F(x) = { F ^ for < x < x s 
K J \ F 1 for x s < x v ; 

and with 

*.= p7^ , F = ^- , F l = 2GM (124) 

in terms of the physical parameters given by the initial radius R s and the total conserved 
mass M; and (ii) quadratic models where 

F(x) = { F °(t-r) fOT0 ^^ (125) 
Fi for x s < x 
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and 



d CM 
x s = s , F = 15— , F 1 = 2GM. (126) 

Xmav ) X~ 



s 



Although we restrict ourselves to these two density profiles, there still are many possible 
sets of values for R s and M. We choose R s = 1 and 0.1 and M = 0.01, where the Planck 
length £p is chosen to be unity. In units used for the numerical analysis, we have a critical 
radius R* = \f^j2l? ~ 0.25£p. So, R s = 1 and 0.1 represent the cases where the initial 
size of the dust cloud is above and below the critical one, respectively. As we have already 
seen, we can recover the general mass scale by rescaling F(x) and t as j3 2 F(x) and t//3. 
In other words, the dynamics of the dust cloud follows this scaling relation. However, it 
should be noted that as the corrected Misner-Sharp mass will not scale in such a simple 
way, the kinematics of null geodesies on the corrected spacetime will not follow this scaling. 
This means that the condition for horizon formation may depend on the mass scale, i.e. 
(3. Due to our choice for the mass parameter M = 0.01, we can see the correction effect 
on the null expansions very clearly. The outer boundary of the calculated region is chosen 
to be -Ro^max) = 2. Our calculation region covers both the classical and effective regimes. 

We have implemented a convergence test to the exact solution for time integration: For 
a = 1 and / = 1, we have the marginally bound Lemaitre-Tolman-Bondi solution 

/ 3 \ 2 ^ 3 

R cx (t,x) = ( i?(0,x) 3 / 2 - -y/F{xjtJ . (127) 

Fig. [2] shows the residual of the numerical solution R n i from the exact solution i? ex (t,x) 

\R n ,i — Rex(t n ,Xi)\ (128) 

for i? ,i — 10 and i? ex (0,Xj) = 10, where n and i label the time step and the spatial grid 
point, respectively. The dust parameters are set to M = 1 and R s = 10. One can see that 
as we decrease the time step At while fixing Ax = 0.05, the residual decreases as (At) 2 . 

In the following we fix Ax = 10~ 4 , where x = 1 corresponds to the outer boundary of 
the calculated region. We have confirmed that the numerical solution will not qualitatively 
change if we double Ax. The time step At is chosen so that the physical quantities on 
each grid point will change their values within 1% at each time step. 



5.2.2 Classical collapse 

Since there is no characteristic scale in classical theory, we have two independent scalings of 
R(t, x), F(x) and t to r]R(t//3, x), (3 2 rfF(x) and t//3, where (3 and r\ are constants. Hence, 
we can recover the results for general radius and mass parameters from a simulation done 
with only one set of parameters. 

Fig. E] shows the collapse of an initially homogeneous ball in classical general relativity. 
The total conserved mass M is set to be 0.01 and the initial radius of the cloud surface R s 
is set to be 0.1. Fig. [3] (a) shows the evolution of the density profile, where the conserved 
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Figure 2: Convergence test for time integration. 
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Figure 3: The collapse of a homogeneous dust ball with R s = 0.1 and M = 0.01 in clas- 
sical general relativity: the snapshots at t = 0, 0.0887925729, 0.125072481, 0.139516828, 
0.145267646 and 0.147557255 are plotted. 
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mass density and the effective density coincide with each other in the classical case. Fig. [3] 
(b) shows the evolution of the velocity profile. Fig. [3] (c) shows the evolution of the mass 
profile, where also the conserved mass and the Misner-Sharp mass coincide with each other 
in this classical case. Fig. [3] (d) shows the evolution of the ratio between the Misner-Sharp 
mass m and the area radius R. This ratio becomes a half at trapping horizons. From Fig. [3] 
(a), we can see that the density profile in the ball remains uniform during the collapse. 
The solution is given by the marginally bound Oppenheimer-Snyder solution [98] where 
the singularity is massive and spacelike. In this exact solution, the singularity appears at 
t = v / 2/3(0.1 3 /0.01) 1 / 2 ~ 0.1490712 - ■ ■. Not only an event horizon but also a trapping 
horizon always appear in this solution. The singularity is always hidden within the event 
horizon as well as the trapping horizon. Cosmic censorship holds in this collapse model. 
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Figure 4: The collapse of an inhomogeneous dust ball with i? s = 0.1 and M = 0.01 in clas- 
sical general relativity: the snapshots at t = 0, 0.0736495616, 0.0898356836, 0.0933232155, 
0.0940745701, and 0.0942363423 are plotted. 



Fig.H]shows the collapse of an initially inhomogeneous ball with R s = 0.1 and M = 0.01 
in classical general relativity. From Fig. H] (a), we can see that the central density grows 
very rapidly while the surrounding region falls into the central region more slowly. This 
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induces strong inhomogeneity near the center and finally the calculation breaks down 
soon after t = 0.0942363423. In fact, the solution is exactly given by the marginally 
bound Lemaitre-Tolman-Bondi solution [2TL l22l [23]. The peculiar behavior at the center 
seen in the numerical solution presents a shell-focusing singularity. In the present class 
of the Lemaitre-Tolman-Bondi solutions, the shell-focusing singularity has been shown 
to be massless, generic and locally naked [55], and moreover curvature strong [58l [59] . 
It can be globally naked depending on the values for R s and M. We can determine 
whether a trapping horizon forms or not before the singularity formation by looking at 
the value of m/R shown in Fig. H] (d), which gives a maximum of about 0.30 achieved at 
t = 0.0942363423. Although this might be slightly larger if we go closer to the singularity, 
the real value is not so different from 0.30. Since this ratio is a half at a trapping horizon, the 
present result means that no trapping horizon is formed in this case before the singularity is 
formed. The existence of a trapping horizon implies an event horizon outside or coinciding 
with it (but not vice versa) in classical general relativity [99J. 

5.2.3 Inverse triad corrections: First version 

Fig. [5] shows the collapse of an initially homogeneous ball in the first version of consistent 
inverse triad corrections from loop quantum gravity with i? s = 1 and M = 0.01. Besides 
the conserved mass density p CO ns, we can naturally define the effective density by 

777- 

= ^gWr>' (129) 

This is defined so that when we integrate this with the invariant 3-dimensional volume 
element on the constant t spacelike hypersurface we recover the Misner-Sharp mass. This is 
directly related to the (t, t)-component of the Einstein curvature tensor and not necessarily 
positive definite. 

Fig. E] (a) shows the evolution of the profiles of both the conserved mass density (dashed 
line) and the effective density (solid line). We can see that the cloud becomes inhomoge- 
neous in spite of its initial homogeneity. From Figs. [5] (a) and (b), we can see that the 
collapse is strongly slowed down in the central region R < 0.05, while the collapse con- 
tinues to take place as in the classical case for the outer region R > 0.2. As a result, as 
we can see in Fig. [5] (a), the conserved mass density at the central region remains almost 
unchanged, while it increases almost homogeneously in the outer region. We can also see 
in the same figure that the effective density is diverging at the center. As the cloud surface 
falls inside R ~ 0.2, the effective density is still nonzero even outside the cloud surface. 
Moreover, for 0.11 < R < 0.24 the effective density becomes negative, which is not shown 
in Fig. [5] (a). A spike develops also in the conserved mass density field at R ~ 0.06 and 
then the calculation breaks down soon after t = 4.56404342. The spike in the conserved 
mass density field should be identified with a curvature singularity, which we can identify 
with the shell-crossing singularity as it also appears in the classical Lemaitre-Tolman-Bondi 
solution. Shell-crossing singularities can be naked but gravitationally weak |100[ll01j . This 
singularity is so weak in curvature strength that it is generally believed to be extendible 
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in a distributional sense J52J 53] . Fig. [5] (c) shows the evolution of the Misner-Sharp mass 
and the conserved mass as a function of R. Although their total values are both 0.01, 
their distributions are quite different for R < 1. The Misner-Sharp mass dominates the 
conserved mass for R < 0.2. It takes a maximum ~ 0.04 at R ~ 0.11. The Misner-Sharp 
mass is a decreasing function of R for 0.11 < R < 0.24, which implies the effective density 
is negative there. Fig. [5] (d) shows that the center is always marginally trapped. This is 
actually seen from the definition of the Misner-Sharp mass m^ 1 ' in Eq. (11031) . Except at 
the center the ratio m/R is less than a half, implying that the shell-crossing singularity is 
not covered by a trapping horizon. 




Figure 5: The collapse of a homogeneous dust ball with R s = 1 and M = 0.01 in the 
first version of loop quantum gravity: the snapshots at t — 0, 3.06904628, 4.04642236, 
4.39347087, 4.51871685, and 4.56404342 are plotted. In (a), the solid and dashed lines 
denote the effective density and the conserved mass density, respectively. In (c), the solid 
and dashed lines denote the Misner-Sharp mass and the conserved mass, respectively. 

Fig. [6] shows the collapse of an initially inhomogeneous ball with R s = 1 and M = 0.01 
in the first version of inverse triad corrections in loop quantum gravity. The qualitative 
properties are the same as in the initially uniform case. Also here, the collapse of the 
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central region R < 0.05 is strongly slowed down, while the outer dust falls onto the slowly 
collapsing central region. Again, a spike develops both in the effective density field and 
the conserved mass density field at R ~ 0.06 and then the calculation breaks down soon 
after t = 3.00299176. The maximum value of m/R is a half attained at the center during 
this simulation as seen in Fig. [6](d). This means that for this case, no trapping horizon is 
formed before the spike or shell-crossing singularity is formed. 

Without loop quantum effects, the collapse generically ends in the formation of a shell- 
focusing singularity. Hence, in regard of singularity formation, the collapse ends in a tamer 
shell-crossing singularity prior to the possible shell-focusing singularity due to the present 
loop quantum effects. 
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Figure 6: The collapse of an inhomogeneous dust ball with R s = 1 and M = 0.01 in the 
first version of loop quantum gravity: the snapshots at t — 0, 2.00754321, 2.65773896, 
2.88917737, 2.97274431, and 3.00299176 are plotted. In (a), the solid and dashed lines 
denote the effective density and the conserved mass density, respectively. In (c), the solid 
and dashed lines denote the Misner-Sharp mass and the conserved mass, respectively. 

If the initial radius of the cloud is smaller than the critical radius R* — 0.25, the 
situation becomes slightly different. Fig. [7] shows the collapse of a homogeneous dust ball 
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with Rg = 0.1 and M = 0.01. As seen in Fig. [7] (a), the conserved mass density has its 
maximum at the cloud surface. Then a spike develops at R ~ 0.046 and the calculation 
breaks down soon after t = 0.156702662. The effective density outside the cloud is positive 
for R < 0.11 but negative for 0.11 < R < 0.24. The Misner-Sharp mass takes its maximum 
value at R ~ 0.11. Outside the cloud, the profile of the Misner-Sharp mass and therefore 
the effective density are almost unchanged compared to the classical behavior as seen in 
Fig. [7|(c). In fact, this feature is also seen for the cases where the initial radius is larger than 
the critical radius. Fig. [7J (d) shows that the maximum value of the ratio m/R greater than 
a half and this is attained at R ~ 0.046, which implies that the shell-crossing singularity 
is covered by a trapping horizon in this case. 
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Figure 7: The collapse of a homogeneous dust ball with R s = 0.1 and M = 0.01 in the 
first version of loop quantum gravity: the snapshots at t — 0, 0.105124094, 0.139369943, 
0.151074789, 0.155215854, and 0.156702662 are plotted. In (a), the solid and dashed lines 
denote the effective density and the conserved mass density, respectively. In (c), the solid 
and dashed lines denote the Misner-Sharp mass and the conserved mass, respectively. 

For the collapse of an initially inhomogeneous ball, which is shown in Fig. [HJ, the quali- 
tative feature is almost the same. However, the central region where the collapse is strongly 



41 



slowed down is much smaller than in the case where the initial radius is larger than the 
critical radius. As a result, a spike develops at the radius R ~ 0.012. The shell-crossing 
singularity is covered by a trapping horizon in this case. 
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Figure 8: The collapse of an inhomogeneous dust ball R s = 0.1 and M = 0.01 in the 
first version of loop quantum gravity: the snapshots at t = 0, 0.134900072, 0.18447396, 
0.202429874, 0.208934676, and 0.211291203 are plotted. In (a), the solid and dashed lines 
denote the effective density and the conserved mass density, respectively. In (c), the solid 
and dashed lines denote the Misner-Sharp mass and the conserved mass, respectively. 



5.2.4 Inverse triad corrections: Second version 

Fig. [9] shows the collapse of an initially homogeneous ball with R s = 1 and M = 0.01 in 
the second version of consistent inverse triad corrections in loop quantum gravity. We can 
easily see that in spite of the very different formulation, the key feature that the collapse of 
the central region is strongly slowed down is the same as in the first version. The radius of 
this almost stopped central region is ~ 0.2, which is much larger than in the first version. 
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Figure 9: The collapse of a homogeneous dust ball with R s = 1 and M = 0.01 in the 
second version of loop quantum gravity: the snapshots at t — 0, 1.43300373, 2.34719999, 
2.95414244, 3.36400543, and 3.39105418 are plotted. In (a), the solid and dashed lines 
denote the effective density and the conserved mass density, respectively. In (c), the solid 
and dashed lines denote the Misner- Sharp mass and the conserved mass, respectively. 
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In the conserved mass density profile a spike develops at R ~ 0.19, t = 3.39105418. 
There is a dip in the profile of the conserved mass density at R ~ 0.25, which corresponds 
to the critical radius R = R*. For R> i?*, the collapse proceeds as it does classically. As 
in the first version, the effective density is diverging at the center. If we go further outside, 
it decreases very rapidly to negative values and then turns to increase to positive values. 
The velocity profile shows two minimum values at the cloud surface and at R ~ R* as seen 
in Fig. [9](b). The Misner-Sharp mass dominates the conserved mass for R < 0.2 and it has 
a maximum about 0.08 at R ~ 0.18 as seen in Fig. [9] (c). The center is always marginally 
trapped and the surrounding region is untrapped. Hence, the shell- crossing singularity is 
not covered by a trapping horizon. 
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Figure 10: The collapse of an inhomogeneous dust ball with R s = 1 and M = 0.01 in the 
second version of loop quantum gravity: the snapshots at t — 0, 0.937882973, 1.5570413, 
2.00350064, 2.2452289, and 2.2463422 are plotted. In (a), the solid and dashed lines denote 
the effective density and the conserved mass density, respectively. In (c), the solid and 
dashed lines denote the Misner-Sharp mass and the conserved mass, respectively. 

Fig. HO shows the evolution of an initially inhomogeneous dust ball with R s = 1 and 
M = 0.01 in the second version of loop quantum gravity. The qualitative features are 
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common with the homogeneous case except that the velocity profile has only one minimum 
during the late stage of collapse. A spike develops in the density fields at R ~ 0.19, 
t = 2.2452289. The spike or shell-crossing singularity is not covered by a trapping horizon. 
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Figure 11: The collapse of a homogeneous dust ball with R s = 0.1 and M = 0.01 in the 
second version formulation of loop quantum gravity: the snapshots at t — 0, 0.550926063, 
1.1729359, 1.83622205, and 2.39522102 are plotted. In (a), the solid and dashed lines 
denote the effective density and the conserved mass density, respectively. In (c), the solid 
and dashed lines denote the Misner-Sharp mass and the conserved mass, respectively. 

Fig. [TT] shows the evolution of an initially homogeneous dust ball with i? s = 0.1 and 
M = 0.01. In this case, the collapse proceeds very slowly in comparison to the classical 
evolution. It should be noted that in classical general relativity, the collapse ends in 
singularity formation at t — y/2/3 x (O.^/O.Ol) 1 / 2 ~ 0.1490712 • ■ • . In the second version 
of loop quantum gravity, the velocity field is kept very small within the whole cloud as 
seen in Fig. [11] (b). The profile of the conserved mass density shows a spike developing 
at the cloud surface at R ~ 0.062 as seen in Fig. [TT] (a). The simulation breaks down 
due to this spike soon after t = 2.39522102. The profile of the effective density is very 
different from that of the conserved mass density. The effective density becomes negative 
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for 0.18 < R < 0.25. This is seen better in Fig. [IT] (c). The maximum value of the 
Misner-Sharp mass is about 0.08. The maximum value of the ratio m/R is a half, which 
is attained at the center. Hence, the spike or shell-crossing singularity is not covered by a 
trapping horizon. 

Finally, Fig. [12] shows the collapse of an initially inhomogeneous dust cloud with R s = 
0.1 and M = 0.01 in the second version of loop quantum gravity. Also in this case, the 
evolution is strongly slowed down. As seen in Fig. [12] a spike develops in the density fields 
inside the cloud at R ~ 0.04, t = 6.4440209. The evolution outside the cloud is identical to 
that in the homogeneous case. The Misner-Sharp mass profile and therefore the effective 
density profile are almost constant in time. The maximum value of the Misner-Sharp 
mass is about 0.08 attained at R ~ 0.18. The maximum value of the ratio m/R is a half 
attained at the center, which implies the spike or shell-crossing singularity is not covered 
by a trapping horizon. 

5.3 Summary of the numerical analysis 

In summary, the numerical results show that the loop quantum effects significantly slow 
down the collapse of the central region which is as large as y/^lp, while in the outer region 
these effects are not important. (There are however quantitative differences between the 
different versions of consistent equations. In the second version the collapse of the central 
region is more strongly slowed down and the radius of this central "core" is larger than 
in the first version.) As a result, the density spike develops at around the radius ^lp, 
which results in the formation of a shell-crossing singularity. Another key feature is that 
the center is always marginally trapped due to the inverse triad correction to the metric. 
This suggests that the shell-focusing singularity appearing at the center might be covered 
by a horizon. 

Since the collapse of an inhomogeneous dust cloud in classical general relativity gener- 
ically ends in shell-focusing naked singularities, this means that although the singularity 
formed in the generic spherical dust collapse in the present formulation of loop quantum 
gravity is still locally naked, the shell-crossing singularity appears prior to the shell-focusing 
singularity and hence the curvature strength of locally naked singularities formed in the 
gravitational collapse is weakened due to the loop quantum effects. In a naive sense, this 
seems to favor the cosmic censorship hypothesis in effective loop quantum gravity because 
shell-crossing singularities are generally believed to be extendible in a distributional sense. 
Moreover, since the shell-crossing singularity will appear near the center, it is likely to be 
trapped due to the inverse triad correction to the metric. We have seen that this is the 
case for some collapse models in the first version. 

However, one should be careful with any definite conclusion because the present work 
can only provide insights into some of the quantum effects in gravitational collapse and 
spacetime singularities. (For instance, as mentioned earlier we have not included holonomy 
effects in the numerical analysis because that would require a more detailed analysis of 
the parameter space.) The slow-down of the collapse indicates that repulsive forces of 
quantum gravity are indeed at work, similar to those that resolve big bang singularities in 
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Figure 12: The collapse of an inhomogeneous dust ball with R s = 0.1 and M = 0.01 in the 
second version formulation of loop quantum gravity: the snapshots at t — 0, 1.50057314, 
3.23870859, 5.12101428, 6.4440209 are plotted. In (a), the solid and dashed lines denote 
the effective density and the conserved mass density, respectively. In (c), the solid and 
dashed lines denote the Misner-Sharp mass and the conserved mass, respectively. 
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homogeneous models. But in our inhomogeneous context this does not appear sufficient to 
provide a uniform bound on energy densities. The more complicated nature of the problem 
is indicated by the formation of shell-crossing singularities which may be a consequence of 
the still present spherical symmetry and the dust idealization used for matter. 

6 Conclusions 

We have analyzed the existence of LTB-type models in the framework of loop quantum 
gravity, starting with an implementation of the corresponding class of spatial metrics at 
the kinematical Hilbert space level. We first discussed the spherically symmetric setting 
and in particular noted the role of lattice refinements in Sec. 13.41 The discussion turned 
out to be much cleaner than in purely homogeneous settings, especially regarding the 
scale dependence of quantum equations. As an immediate consequence, not all refinement 
models used so far for anisotropic homogeneous models can be embedded in spherically 
symmetric models. In particular, non-trivial refinement schemes which would give rise to 
equidistant difference equations (where point holonomies are of the form exp(i/i(p / )c/) with 
only the variable p 1 conjugate to a connection component cj entering) do not seem realizable 
in spherical symmetry. It is thus important to improve the analysis of non-equidistant 
difference equations, for instance along the lines of |1021 1103] . Despite some restrictions on 
the refinement scheme, a whole class of refinements varying with the spatial size remains 
allowed. All these schemes have the correct scaling behavior under coordinate changes, 
even though further reduction to homogeneous models may suggest improper scalings for 
some of these models. This shows that it is only the reduction to homogeneity which may 
suggest improper scalings because scaling- dependent parameters arising in the reduction 
are overlooked. The restriction of refinement models based solely on their rescaling behavior 
in homogeneous models is thus too strict. 

We then turned to the LTB conditions and showed that the classical conditions translate 
easily to the kinematical quantum level, which allows further studies of quantum reduction 
mechanisms of loop quantum gravity along the lines of [7] . However, the constraint algebra 
makes a discussion difficult at the dynamical level where the consistency issue of the op- 
erator algebra of constraints together with LTB conditions would have to be analyzed. At 
this stage, further progress is possible based on an effective treatment of correction terms. 
While we have not derived complete effective equations which would contain all relevant 
correction terms at once, the inclusion of individual correction terms of certain types can 
already be used to see how quantum effects can be realized consistently. 

Indeed, we have found consistent formulations of the LTB reduction for different types of 
corrections: the constraint equations (|66|) . (!73|) and (|92|) with the evolution equations (1671) . 
(!74|) and (1911) . respectively. Moreover, the effects can be combined in consistent equations 
summarized in Sec. 14.21 The consistency of an anomaly-free formulation relates corrections 
in different terms of the Hamiltonian to each other, but also to required corrections of the 
classical LTB conditions. Here, our general understanding of an LTB-type reduction is that 
metric coefficients in ([2]) are related to each other by a proportionality L oc R' . Classically, 
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the factor is one, but consistent quantum corrections require a non-trivial dependence on 
R or R which would be important in some regimes. Thus, while we have the same type of 
reduction of degrees of freedom, explicit dynamical implications for the metric may change. 
In particular, the relation between L and R affects, for instance, the position and behavior 
of horizons in addition to what a change in the dynamical behavior of R would imply. Also 
the behavior near classical singularities changes due to correction factors in the metric 
which may even diverge right at a classical singularity. A complete space-time analysis of 
the resulting effective metrics remains to be done. 

The change in the form of the metric had unexpected implications: Further reduc- 
tions which are possible classically, such as Friedmann-Robertson-Walker solutions, no 
longer exist, although on large scales there are approximate such solutions. This is quite 
unexpected and suggests that caution is necessary regarding the dynamical realization of 
homogeneous models. Moreover, in combination with corrections to the classical equations 
of motion, we have seen new terms in effective densities which can become negative even 
where mass densities remain positive. Also numerical simulations suggest that repulsive 
forces of quantum gravity are active on small scales. All this can affect the horizon behav- 
ior as well as singularities if negative energies and corresponding repulsive forces become 
strong enough. However, our analysis of central singularities, which may be spacelike 
or non-spacelike, did not reveal any indication that they would be prevented completely. 
While spherically symmetric loop quantum gravity is singularity-free at a fundamental level 
of difference equations [T2] (see also [1041 1105] for Gowdy models), the development of in- 
tuitive geometrical pictures requires non-singular equations for an effective geometry. We 
have not analyzed full effective equations and did not yet consider all possible correction 
terms; it may be that some of the corrections which are more complicated to include can 
avoid singularities more generally. Nevertheless, the difficulties in avoiding inhomogeneous 
singularities, in contrast to the fact that several spacelike singularities have been shown to 
be prevented based on homogeneous models of loop quantum gravity, suggest that general 
singularities present a qualitatively different issue compared to what has been realized so 
far. From our numerical simulations we might speculate a possible resolution mechanism 
more subtle than the analog of a cosmological bounce: Strong curvature singularities seem 
to be replaced by weak singularities which may be extendable by distributional solutions 
or with more realistic matter models. 
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